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Q» 1 ABSTRACT 

The general equations governing the large amplitude 
flexural vibration of any thin elastic shell using curvi- 
linear orthogonal coordinates arc derived in this report. 

These equations consist of two coupled nonlinear partial 
differential equations in the normal displacement w and the 
stress function F, From these equations, the governing equa- 
tions for the case of shells of revolution or flat plates can 
be readily obtained as special cases. The material of the 
shell or plate is isotropic and homogeneous and Hooke’s liw 
for the two dimensional case is valid. Modal equations are 
derived on the basis of a series for the assumed vibration 
mode. A typical term of the series is taken to be separable 
in the space coordinates and time. The nature of the modal 
equation obtained for the case of flat plates exDlains the 
reason v.’hy the nonlinearity associated with the large ampli- 
tude flexural vibration of thin flat plates of any geometry 
is of the hardening type. Similarly, the modal equation 
applicable for the case of any thin shell shows the reasons 
why the softening type of nonlinearity is quite possible in 
the case of thin elastic shells. The difference betv/een the 
hardening type of nonlinearity in the case of flat plates 
and straight beams and the softening type of nonlinearity 
in the case of shells and rings can in general be traced to 
the curvature (or lack of it) of the undeformed median surface 
of the structure concerned. 
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Introductio n. 

The nonlinear (large amplitude) flexural vitrat-'n 
of plates and shollo have received considerable attention in 
the last 20 years. This area is of great importance and 
interest in the case of structures of low flexural rigid i;y 
which as a consequence are easily deformable. The import; oicc 
of this investigation arises from the fact that when the 
flexural vibrations Involve large amplitudes, the frequency 
of free or forced vibration is very much dependant upon the 
amplitude. Significant changes in the frequency correspond- 
ing to the small -amplitude vibration can occur when tho 

amplitudes involved aro no longer small so that the resulting 

« 

governing equations of the problem are nonlinoar. 

A somewhat comprehensive bibliography of recent 
papers published in tho field of large amplitude flexural 
vibration of thin elastic flat plates and shells is given 
at the end of this report. This bibliography which gives 
82 references on plates and 57 references on shells is taXen 
from a recent survey paper (to be published) titled ’Nonlinear 
Flexural Vibrations of Certain Deformable Bodies* co-authored 
by this author (Ref. 6), 

As can be seen from the bibliography mentioned 
earlier, the existing literature as far as the nonlinear 


flexural vibration of thin flat plates is concerned, relate 
to plilco of rectangular, circular, elliptical, triangu^* 
and parr.llelograirjuic geomotry. Moot of the investigations 
carried out so far deal with inotropic platoo. In order to 
obtain the solutions of problems of plate geometries other 
than those already investigated, the governing equations for 
the large amplitude flexural vibration of thin elaotio flat 
plates in orthogonal curvilinear co-ordinates are necessary. 
Therefore these equations have been derived in the next 
section where the necessary modification to tackle the problem 
of buckling of thin clastic flat plates is alco indicated. 

In the area of shells, the existing literature 0:1 the 
large amplitude flexural vibrations relate to circular cylin- 
drical shells, spherical and conical shells. Recently under 
the guidance of the author, work was carried out on the non- 
linear flexural vibrations of orthotropic oval cylindrical 
shells. It is of interest to study the behaviour of thin 
elastic shells of geometry other than those mentioned above. 

To this ond, the general equations governing tho nonlinear 
flexural vibration of thin shells of revolution have beon 
obtained in Section 3 . It is proposed to use these equations 
in the case of the ogival shell., torus shell and ellipsoidal 
shell of revolution. The governing equations mentioned abov* 
are suitably modified to tackle the problem of buckling of 
thin elastic shells of revolution. 


Having derived the governing equations for the non- 
linear ilexural vibration of shells of revolution as well aa 
of thin elastic flat plate3 using curvilinear orthogonal co- 
ordinate a, the governing equationo for tho nonlinear flexural 
vibration of ANY thin elautic Bhell can be readily obtain ;d. 
From these very general equations-, the governing equation 3 
applicable for ohells of revolution or thin flat plateo cun 
be obtained as special cases* 

In the last but one seel ion, on tho basic of on 
assumed vibration mode of the product type (product of func- 
tion of time and function of space variables), the modal equa- 
tion is determined for the genoral case of thin clastic shells. 
The modal equation is then specialised for tho case of flat 
plates end beams. The characteristic features of the modal 
equation as applicable to thin shells on tho one hand and thin 
flat plates and straight beams on the other hand are discussed. 
It is conclusively shown that in the cace of flat plates .ind 
beams, the nature of the modal equation is such that the non- 
linearity of the hardening type results. 

In the last section, the governing equations for the 
postbuchling behavior of flat plates and thin shells obtained 
as a byproduct of the study of the problem of large amplitude 
flexural vibration are examined. It is shown that in general 
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the set of algebraic equationa in tho coefficicnto of the 
assumed series for the postbuc 1-ling mode arc a oet of 
coupled cubic equations in th'-se coefficients in the caoe 
of both platea and shells. However, in the absence of any 
lateral load acting on the flat plate, the set of algebraic 
equations in the coefficients reduces to quadratic equations. 
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In thio section, the governing equations for tho 
large amplitude flexural vibration cf thin elastic flat platen 
in curvilinear orthogonal oo-ordinalea 0c : and cx 2 are 
derived. , <*2 and ^ form as uaial a right-handed ortho- 
gonal coordinate system. 

An elemental length do on tho flat plate median 
surface is givon by 

(ds) 2 = (A, do6j) 2 + (A 2 da 2 ) 2 (2.1) 

where, end A^ are Lame^ parameters of the median surface aid 
are functions of the curvilinear orthogonal coordinates cL x aid a 2 . 
For a point off the median surface at a distance ^ away, the 
extensional and shear strains are given by [l] , [ 2 ] 


€11 = ^1 + ^ ^11 

€ 22 « € 2 + ^X 2 z 
V*. = « +2'3fX l2 

( 2 . 2 ) 


where, the median surface strains 


€,, e 2 and 05 are given by 


t 6 i 





03 



aco/A 2 ) 

<3oCi 


Aa (3WAi) 
A'J 3a. 2 



(2.3) 


The curvature expressions X^X^and X 1Z arc given by: 


^11 - 
'Kn - 


J_ a_ A_ 8w \ 1 /3A, \/l_ 9»V \ 

Ai 3(tj \Aj detj/ Aj A 2 \dtt-2' vv2 3ot 2 / 

9 

_1_ _a_ /^l_ 3wV _A_ ^AaVl 3 k/ \ 
A2‘ d a 2 V^2 3oC/j/ A,A2V3 oc|/U 1 doc/ 


v [a 2 * _ jl 1 _ 1 /a^yawV 

• A,A2 3w ( 3a 2 A t \3«j/\3ot|/ A2 \3fc,/ V3«-2/_ 


(2.4) 


u, v and w are as usual the components of the displacement 
vector in the orthogonal co-ordinate directionc ot, , orj. and 
respectively. Prom equations (2.4) it follows readily that. 


-V 2 W = 



(2.5) 
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For the isotropic case, the strain erercy u can ho 
written an (See Equation (9*12) cf Ilef.2): 

u = u, + u 2 (2.6) 


where , 


ll, - Extensional otrain energy 


l tiV -2 0-V)(e,e 2 - “-)]dA 


(2.7) 


and 


u 2 - Bending strain energy 

- -£* ^{^i| "t ^ 22 ) “ 2 C*- V) [*„ ?'22*"(^/2^ ]] 


( 2 . 8 ) 


In the above, C = 


Eh 


and D = 


Eh 3 


0-V 2 ) 12 ( 1 - V 2 ) 

In view of Equation (2.5), Equation (2.8) can he written as 


U 2 = | J A {(7 2 W) 2 -2 0 -y)[X„ 7C n -C«| 2 ) 2 ]}dA 


(2.9) 


In the above integrals, the elemental area dA of the 
median surface of the flat plate is given by 


dA ** A, A2 d^, d « 2 


( 2 . 10 ) 


0 I 


From the principle of minimum total potential, &(u-VJ)u: ) 
where the work done by the applied lo ode is given by 


W = J A ( ( i >u + (3° +Zw )d.\ 


( 2.11 ) 


In the above, j>,S) and Z are the components of the applied 
load in the coordinate directions • od, , of 2 and ^ respec- 
tively. The three equations of equilibrium obtained from rhe 
• principle of minimum total potential can bo shorm to be: 


_a_ 

da, 


(A 2 Nii) - N22 (f^) + aT 2 n (2 ) •» A, a 2 <£ =0 


(2.12) 


§k( A ' N 21 ) - n„(|^V (A| Nil) *■ A, A 2 © = 0 


D[(»,.0 V‘»- ( ,-V>l] g] 


+ 2.N 


(2.13) 


'©g '©£}-*'*** ■ » 


In the above, V*w * v 2 (v a w) and M^, N and N 


22 


12 


<4 

are the stress resultants defined by N- =f h Oij d%. 


A3.su, the. function L(A,,A 2 ,w') stands for 
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L = 



\ 

(2.15) 


where, the curvatures and K n are defined by 

Equations (2.4) in terms of the displacement component w. 


It is of importance to recognise that Equations 
(2.12) Jind (2.13) are the equations of inplane equilibriua 
and are exactly the same as those developed by Tiushtori In 
the well-known Mu ditar 1-Vlacov approximation of thin shell 
theory (See Section 17, p.84 of Ref .12)). 

The expression inside the brace *(•••• } appear- 
ing in Equation (2.14) can be simplified by means of Equa- 
tions (2.12) and (2.13) wherein the load components and© 
are neglected. Thus Equation (2.14) can be written as: 


~ L DO-V) 
DV 4 W- — — 
A, ^2 


.m rj. ±jl hi\j±-\±hi\(± 

11 A| 5<X|\A, ''A,A2^3 a 2 A 2 9^2'. 


-nJI - a - (±- — Y] 

z LAz B*2 \A 2 Baa' U,A/\ ax, M, d«yj 

“2N.*[ 


dcC| 0)^2 


/ Sa* y 

\ dvv \ 

\ a<c, A a 

kj ^<2 - 
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In view of Equations (2.4), the above equation can also bo 
v/ritten to 


dv 4 w- 


DO-V) 
A, A 2 


L + N|| 96, ( t N 22 9^2 + 2 N 12 ?h 2 



(2.17) 


The function L given by Equations (2.15) if 
written out purely in terras of W will be somewhat conipli- 
* cated and lengthy. However, under certain conditions - to be 
discussed below - this function .‘0 can be shown to be zero. 
For example, if the Lame* parameters A^ and A^ are such that 

the products of the derivatives of A^(or A^) with 

7 - ~ (or 7 - ) or their derivatives con bs 

A, aoc, \ A z doL Z J 


neglected in comparison with the terms of the type 

±.(L **-) or — (— — ) 

3oC|\A| d^i/ dotj \A2 dcX. 2 ./ 

then, the function L defined by Equation (2.15) can bo 
shown to vanish. In such a case, the curvature expressions 
given by Equation (2.4) also must be simplified to read: 
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1 d /I 

3w\ 

v 1 

3 / 

<\ 3w \ 

Aj 3<x f \ A, 

IZj ’ 

* u ^ "a. 

doc^' 

^2 dtf.2/ 

%I2 

1 

d 2 W 



1 *'* 

l 





( 2 . 18 ) 
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Consequently Equation 

(2.17) 

simplifies to 


DV7%, - N "U- ki 

1 Y 

_N “[i iSi 

• --T1 

A, 3ot,/. 



d*W ] 

n 

1 



(2.19) 


Another way to look at this approximation (that is, L = 0) 
is to consider the expression for the strain energy U. a 3 
given by Equation (2«6) from which it is seen that U. con- 
sists of two parts, one proportional to the extensional 

Eh 


rigidity C «* 


(l-v 1 ) 


and thj second part proportional 
Eh 3 


to the bending rigidity L> =* • Prom Equation 

(2.9) it is clear that the bonding strain energy U* 


2 . 2 - 

depends upon two quantities, namely, (V Z W) 


and 


2. O V> [»£„ Xn-Qin.) 1 -] 


If the integral of the 


term 2 C l- V) [x „ X zz - C K /2 ) 2 J is neglected in 

comparison to the integral of (V 2 W) 2 , then d A 


and in such a case also, L will vanieh. 

The stress rcsultantf N. IT and N are func- 

11 22 12 

tions of the stress function F( oc, , <* 2 . ) as defined later 
on so that Equation (2.19) represents an equation involving 


12 i 


W and F. If the terms involving coupling between V/ nnc: F 
are neglected, Equation (2.19) reduces to the well-known 
equatior. of the linear bending theory of isotropic flat plates, 
namely, PV^w^z • 

It has already been mentioned that Equations (:i.12) 
and ( 2 . 15 ) are the equations of inplane equilibrium and that 
they arc the same as those developed by Muohtari in the 
Mushtarl-Vlasoy approximation of thin shell theory [ 2 ] 

* The stress function F( <x, , ot 2 ) is so defined as to satisfy 
Equations (2.12) and (2.13) exactly and the relationship between 
the stress resultants (N , N jurxd N ) and the stress fimc- 

n 22 12 

tion F are exactly the same as those defined by Vlasov ;.n 

/ 

the Mushtari-Vlaoov approximation of thin shell theory (See 
Equatior. (17.S) of Ref. (2).). Thun, 


N _ _ J_ _d_ f\_ dF \ I /dAA / I £\ 

11 A2 ^0c2 doc 2 / AjArj, \ doc, A A| doc,/ 

N U n - -L — f— —V — f^Y- —) 

A, da, V A, da,/ A,A 2 \da-2/'A2 d<*2/ 

N ”• J— V^yi aPW^Y 1 3 FV] 

12 A,A 2 Wo^Aa, do:,/ \doc.AA 2 doYJ 

( 2 . 20 ) 

obviously, 


(N„ >N 22 ) 



( 2 . 21 ) 
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Because of Equations (2.20), Equation (2.16) or 
(2. IS) represents one equation involving coupling between 
W and ?, Before we proceed to obtain the eocond equation 
involving \7 and F, it is worth noting that Equation (2. IS) 
represents the governing equation of plate buckling provided 
the signs in front of N 22 ani are changed from minus 

to plus in Equation (2.19). Thin in obviously because for 
buckling pr obi ens, and N 22 will be compressive (and 

similarly for N^) 60 that Tor plrfce buckling, the governing 
equation becomes: 


. .1 I £ /I dW\l I 3 /I 3W\" 

DV 4 W + NJm ) +hJo9 ■" — — ( — I 

L A, \/\\ 300,/J LA 2 . 3ot2 VA 2 . 5 0 ( 2 ./. 


-t 2 


n 1 = 7 

,2 La,a 2 a«,da2- 


(2.22) 


If a more accurate equation than Equation (2.22) io required 
for the 1 late buckling problem, we can use Equation (2.16) or 
(2.17) wherein we change from mini 3 to plus the oigns in front 
of the terms involving N^, N 22 ar.d In these equations, 

N^, N go and li 12 represent the stress resultants (acting on 

a typical plato cloment) induced ty the compressive edge loads 

avvd 

that give rise to buckling ^are del errcined from the inplano 
equilibria;!! equations, namely, Equations (2.12) and (2.13). 
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(. 



To go back to the* large amplitude flexural vibrat .on 
problem of the plate, besides Equation (2.16) or (2.19), 
another equation involving \7 mid F has to be obtained. ;‘or 
this, the compatibility equation has to bo brought into the 
picture. The required compatibility equation can be obtaii cd 
by specialising the compatibility equation of Mushtnri-Vlcu ov 
approximation of thin shell theory (See the third of Equations 
(5.1) of Ref «(2) .). % The specialisation required is to let tin. 
radii of curvatures r 1 and r tJ of the shell median surface tend 
to infinity so that the shell medi.in surface flattens out end 
becomes the median surface of the flat plate. Thus, the com- 
patibility equation can bo shown to bo: 




£ _ Aj dl£ _ c)A| j I 

1 dot, 2 doc-j m do^J/ 

e 2 _ k. *£ _ w ¥l}\ 

dd2. 2 d*Cj doCjJJ 


= 0 


12 . 25 ) 


From the stress-strain la* for the icotropic plate 
material and the definition of the stress resultant, it readily 
follows that 


c CN lt -V N zl ) ; C 2 


(2.24 ) 


where the stress resultants , N ?0 and N 12 arc defined in 
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terns of F by Equation (P.,20). Upon introducing th* above 
relations into the left hand side ox the compatibility equa- 
tion (2.23), the resulting enpren.iion after algebraic simpLifi-. 

cation ard manipulation con bo shown to be A 1^P, , A-,\ v/iere, 

Eh ' v x ) 

V (“’ ) ** V 2 [ V 2 (-* )] and the Luplocian operator V 2 

# 

is defined in Equation (2.5) or (:?.2l). In fact, the expression 

A 1 A p 4 

““ErT* ^ V 3?) can also be obtained by specialising (i.c. lot 

* fnd r^ -m) the final eqjation obtained ffcom the com- 

patibility equation for the case if Muohtari-Vlasov approxlma- 
tion of thin shell theory. (Sco S action 17, p. 84-90 of Koff2).). 

The next step is to intr educe 6 f G^ and 60 as given 

by Equation (2.5) into the left hind side of the compatibility 
equation (2.23). In doing this, it helps to remember that 
e V 6 2 as given by Equation (2.3) have a linear part 

and a nonlinear part. If for the sake of convenience and 
simplicity, the nonlinear part is neglected for the time being 
and the resulting linearised relationship for G^, 0^ and o) 
are introduced into the left hand side of Equation (2.23), 
it \7 ill be seen that Equation (2,23) is automatically satis- 
fied. In checking this, use has to be made of the following 
condition of Gauss specialised for the case of the median 
surface if the flat plate by letting the radii of curvatures 
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and of the shell median surface tend to infinity; 


J>_ n 

■i V A i 


3it, l A, dli.J 


d r\ 3A,\ 
3« 2 \A 2 ha 2 / 


(2.2‘j) 


If only the nonlinear ter*ao in the strain-displace- 
ment relations Given by Equation (2*3) are retained and they 
are introduced into the left hand side of Equation (2.22), 
the resulting expression will be: 


_£ f — f Al (— ^L\ 2 l - -L ^V1 aw V 

3tx, LA| Laoc,i 2 VA 2 3aJ J 2 \duc,/ VA, doc.,/ 


i /awyawyaAAil a 1 f j__ /awyawvi 

A, A 2 \3o£, ' V 3o( 2 A 3oC 2'- •> aoc, d« 2 . A,A 2 '3oc.,' Vda^J 


a / j_ r a_ f a, /y aw\ 2 1 _ j/dA, vj_ ay 

doc 2 L A 2 l 3 oc 2 l 2 ^A, doc., / J 2 3oc2''A2 3ot2./ 


A, A 


l M 2 


A3wyawY3A 2 v]i 

Vdcc , / V3 cc 2 / \da , )\ i 


(2.26) 


Ao mentioned earlier, the left hand side of the ccm- 


patibility equation (2.23) is also equal to Al /? / v^F). Jcnce, 

Eh 


A. A 


^4 p 4 

( V F) a Expression given by (2.26) above. (2,27) 


s 17 s 


For working out cry practical pro *c 1cm, expression given by 
(2.26) neeuw to be simplified. Tc thin end, we employ the 
same ap^r jxirration that in involved in making the fimeti.u 
L appearing in Equation (2.H) zero, namely, neglect terms 

that involve product. c with or 

or ot (^2) * the expression 

given by (2.26) simplifies to: 

_J f / 8w \ 2 _ /3 2 w \ / a* tv \| 

A1A2 ia< 1 c£ ( ac^/ Vdoc ( z / i a«f- /„ 

fhu3, Equation (2.27) reduces to: 


V 4 F = Eh 




d 2 W \ 2 


AjA£ a^tj 80(2. 


\ 2 / 1 s z w \/ 1 a a w\) 

/ \a* aocf Aa| aaiJ] 


(2.28) 


We can arrive at the right hand side of Equation (2.28) 
in another way, starting out with a simplified version of the 
compatibility Equation (2.25) wherein the terms 


*“ ("-:&) 
0 


n re neglected in comparison with the terms 
’^c”) similarly, the terms 


and 


A 


doc 2 ) 
with the terms 


a\'e neglected in coL.par5con 


(-S 9 

r*. a *- 


The 


basis of such a simplification is the same as before, namely, 
the termn involving the product of the gradients of tho Lam® 
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parameter. r. and are noglectc d in ccnpcri.von with tc.v '« 
involving the products of the Zone parameters themsolveu* 
Thus tnc compatibility Equation (' .22) simplifies to: 


A 2 d%. A, oC, 

— 1 1 f — - ■ — . -i- * — ■-■■ ■ 

A, dttf A 2 <*« 2 dec, dec 2 


= 0 


(2.29) 


Introducing into the above equation, the nonlinear part of the 
strain-dl o placement relations, namely* 6, - ~ ^t*) 2 ? 

* u ~ 1 (1 Jw\ 2 

2 2 \Ao doc J 


and 


CO = • (*H 2* \ 

A,A2\3Ct, ddiJ 


and neglecting tcrina involving the products of >(^r) 

etc., we will get the right hand tide of Equation (2.23). 


Equation (2.28) representt the most simplified ver lion 
of the compatibility equation and together v/ith the equilibrium 
Equation (2.19) constitute a pair of coupled, nonlinear partial 
differential equations in the two functions V/ and F« The ic 
two equations as they stand arc vt lid for the static caoa. 

For the vibration problem, the necessary modification based on 
the D Alembert's principle has to be mmle. Thus, the loading 
term Z appearing in Equation (2.19) has to bo replaced by 

whero P 10 tho BU»I 

per unit area of the plate and t is time. For free vibration, 
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It is worthwhile suaLiarii ing the governing equations 
for free undamped nonlinear flexural vibration of a thin 
elastic flat pl&te. These equations in curvilinear ortho- 
gonal coordinates a ^ and « 2 arc:- 


V 4 F = Ell 


U 


3 2 w \ 2 yj_ 3 *vy\/ i 

A 2 3ot, \A* daf/\A| d&| /) 


(compatibility equation) (2.30) 


P» + Dv“w4^ 

L A| Doc f \ A, Doc | / A z 'Ai Datj / A ( A Z dat,3ocj. J 


(Equilibrium equation) (2,31 ) 

where, Iho stress resultants ond H are define! 

in terms of the stress function P by Equation (2.20), 

However, in view of the approxlw ition involved in Equations 
(2.30) »ind (2.31 )i these stress resultants can be approxi- 
mately expressed by 



-L S?L 

A\ dec* ’ 



J_ d«F 

a* act* 


N.2 * — 


i a*F 


AjA 2 doCjictj 


( 2 , 32 ) 


The above approximate expressions bring to mind the defini- 
tion of Airy Stress function of two dimensional elasticity. 

A 4 ? 0 

The bihinoonic operator v ic defined by v » v . v 
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2 

where, tho leplocian operator v is given by 

v 2 c • > = -L- (4-rr^ ^1 + \t±) sti 

A, A 2 L c5cc, L\ A, / 9ct, J ScxjIAAj/ Sotj 



(2.33) 


'x'ho above equations can bo very readily specialised 
for the ease of any orthogonal system of coordinates a end 

a 2 » example, cartesian, polar, elliptical co-ordinatee etc. 


With ds, * A, dot, and c s^ = Aidoti 
Equations (2.30) and (2.31) can b-j written as, 


and 


p(v +VC7*W) - {n„ ^ ♦ N* ♦ 2 N« 


where y 


V 2 .^ and 


v 2 C") -/f 1 

lAj 5s, L 9s, J A, ds^L 952 J / 

r a 2 c-) 9 a (-)~ [ 

~ L asf + as| J 


(2.34) 


(2.35) 


2.36) 


3s| 

0Ao , 3 A . 

since the terns involving the proiuct.s of — £ and L 


da 


1 


t)a< 




• PI • 

• tm I » 

arc neglected in the approximation employed here. 

With suitable modifies dons, the pair of equation 

(2.30) ar.d (2.31) [ or, Equations (2.34) end (2.55)] tiro 
represent the governing equations of the problem of post- 
buckling behavior of thin elastic flat plates. For this, 

it has to be recalled that both \7 and F ore independent of 
time t so that the inertia tern ' pw «= 0 ; also, the signs 

in front of the terms involving toe stress resultants H , 

« 

Mg 2 defined in terms of tie stress funotion F by 

Equations (2.20) have to be chong3d to minus. 

Before concluding this section, a brief discussion 
of the methods of solution involving Equations (2.30) and 

(2.31 ) in in order. One method is to assume a deflection 
function \l that satisfies all the boundary conditions of 
the prob’.em. Based on the assume 3 W, the compatibility 
Equation (2.30) is used to determine the corresponding stress 
function F. In general, the assumed W and the correspond- 
ing F will not satisfy oxactly the equilibrium Equation 
(2.31 )• An approximate solution of this equation involves 
using th $ Galcrkln'u method. The details of this method os 
applied to nonlinear flexural vil cation problems are given 

in [3] and [ 4 ] . 


: 22 : 


An alternate method of solution is bar^d on v.- itir-r 
the in plane equilibrium Equations (2*12) and (•-♦• ’‘j ) in Uo ■ : 
of the displacement components u, v and w. For this, :.t 
will be recalled that N„ (£,+ 1 ^ 2 ), n;C| '^ 

and “ i2 = r ^co t v:h »re , 6 1 f e g and co are G ivcn 

by Equations (2*3) • Thus, Equation (2.12) can be written as: 

/ JL |7A*\ ay _ _i_ /9Ar^ a u + _!_ 

latt, L\A,/ du-,J A,A 2 \Sot,> 2 A, SfiCa.l.VAj/daaVA,/! 

♦£[*(£)»] 

ran /gA>v » /3 a ? ,\c>u + a 2 *> 

+ ^1 aoe.L A.laocJ aAsoc./cXX, doc.doc* 


ra fi /aAa\ _ i /a M c’u a^_ 

+ ^l doc,L A, Vaoc,y A,Vaoc,/ £<a., doc, 

a r/A* \ a /u.\| v ^A 1 yaA 2 > \ 

" 2A, a«2i-VA 2 /a«2\A,/J A^vaocxAdct,/ 


-~L -L[a,a2 — 

2 A, a* 2 L VA2 ^jJ 
2 \ 3 <x/ 'A 2 aoc 2 / a&ilA 2 / 'Aj aoL,/ J 

1 9 [A, /awvBwVI JA; /_!_ 3W j 

2A ( 5« 2 l- A 2 'SttjAaflLj/J 12 vA 2 act 2 / 


i pA 2 y i aw f I a f(M 
" ZW.AA, aoc.y 2A, aoc 2 KA 2 y 

w jS^ Page B 


r3W\/ &W V 
^dOC/Vaot^. 


(2.37) 
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In a similar manner , Equation (2, 4 3) can aleo Ki written 
in terma of u, v ar id w# However, it ie easier to obtain 
tli at equation from Equation (2*37) by interchanging 
and A 2 , c:^ and a and interchanging u and v. Having 

written the two inplan ft equilibria equations in terma of 
u, v and w, the next step io to determine from theue tv/o 
equationn, the u end v functiono for on assumed \T* The 
details of thia method of approacli are given in [ 5 ] • 
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Section *»> 

In thie section, the governing equations fox* thy 
large amplitude flexural vibration of thin clastic shells of 
revolution are derived. There equations can thon bo specia- 
lised for circular cylindrical, conical, spherical, ogival, 
torus shells etc., which are all cneciol cases of the shell 
of revolution. 



Figure I shows the meridional curvo which when 
rotated about the axis of revolution generates the shell o:‘ 
revolution. The two principal radii of curvatures of the 


shall median surface are r 1 raid :c 2 and their directions :vb 

any point of the shell median curfr.ce are ao&uncd to coincide 
with the diroctlono of the orthogonal co-ordinates (} an d C. 
For n shull of revolution, ani r 0 ere functions of (j only 

and not of e. Also, no can he seen iron Figure I, 


r 6 tz r z sin f 


( 3 . 1 ) 


and, 


given by 


d h c 
d <}> 


r, cos^ 


The material lav/ for the orthotropio case is 


* 7 i = E|| 6|, + E12 622 

^22 = E -22 ^2Z + ^ 12^11 
r /2 = 






( 3 . 2 ) 


( 3 . 3 ) 


where E |1t E^, E^ and G are four independant materiel 
constants of the orthotropic material. For the isotropic 


case 


E 


II 


E 


E 


22 (1 -V 2 )’ 


G = 


20+D) 


and £,0 = 


EP 


12 (\-v z ) 


The strain at a general point (Cj, 0, ^ ) can be 
written in terms of the 3ix median surface deformntional 
quantiti es e 2 » co , ^ , x £2 and # 12 ac follows: 

e n - 61 + 


^2Z~ 

y, 2 - I0i*2/^X, 2 


( 3 . 4 ) 
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t • 


v/hex'c, the (| - 8 - 2r co-ordinates form as usual a r1c.ht 
handed co-ordinate eye ten. The 'yr co-crdlnnto is pocitivo 
in the direction of the unit n'vi.-u'&l outward an in [sj * 

The median surface deformation cliaractcrioed by the exten- 
eional strains 6^ and and the eh ear etrain. (0 e.a well 

ao the curvatures and X arc defined in 

11 .22 12 

terms of the displacement components u, v and w (which are 
in the (), 0 and y co-ordinate tirectiona reepeotively ) ns 
follows [l] , [ 2 ] 


€| = i- a >+ + ? 


62 = r/»» + \ Cot * + 'r* + i(t w,e ) 


CO = 


P6„ - 


Lv 'a ^ r Lt . , cw,fl)(w,j) 

t) ^ »(, a,s r 2 C fc * :rT 


- 1 


r i »o 


r , a ^ u , 


^22= "7F W ^09 “ 


— (W,*) Cot ^ 
r /'2 


X 


- _ JL i_ /Ll 


12 1: a* 


fi w ) s iL )C!i£.i w 1 
\r 0 V r 0 [_' w ’ 6 ' r„ r, ’<MU 


( 5 . 5 ) 


It is worth noting thrit the exprocoions for 6^ t and a> 
are nonlinear whereas those for X )0 ^ ^ are 


linear and dependant only upon ft and its derivatives, 
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but not on u or v* It in easily seen from Equati one (5.5) 

that 


(*,, * X«) = -V : 


w 


= --{ 

r,r„ l 


Mf-Hi *©"•») 


(3.6) 


The expressions for stnin energy U, extensions], 
strain energy and tho bending -strain energy U 2 are tbe 

Game as those given by Equations (2.6), (2.7), (2.0) and 
(2.9) v.hsrc, the elemental area 1A in tho case of the shell 
of revolution is given by 

dA ■ r 0 r, d^ do 

/. 

The expression for worlc done by the applied surface force 
components $ * ® and Z (per imit area) in the co- 
ordinate directions (}, 0 and nj respectively is 

W = J Atea (^ >u - + ® l/ + 2w) dA 

Using the principle of minimum total potential, namely, 

6 Cu-w) =■ 6U TOt = 0 , the three equations of equili- 

brium ii the Cj , 0 and ^ co-ordmato directions are obtained 
in terms of the stress resultant*. 11^ and the moment 
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reoultnntn which are defined no follows: 

r h /2 


= 4^ 

-h/2 

M ij “kvoijdlf 


(3.8) 


Thus, 


n,| - C (6, •ti’e*') ; n 2 2 - C (c 2+Ve,); (5.9) 


and, 


M„ = D(j£„+?'X2i); M 22 = p(?t2 2 + M,,*D(l-)»K,a (3i10) 


where , 


C = 


Eh 


(\-y 2 ) 


and D = 


Eh 3 


120 - 7 *) 


is the extennional rigidity of the shell 
is the bending rigidity of the shell. 


The three equations of equilibrium are: 


hk° N ") + Je( r ' Ni2 )“ ( r ‘ Cos 0 ) Nn + r ° r > $ = 0 


(3.11 ) 


A(r, n 22 ) + ~-(r„ N e J t(r, Cos (ij Nj 2 + r„ f, 0 = 0 


(3.12) 
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(3.13) 


In tern:! of the stress and moment resultants, the last 
equation can also be written as: 



(3.14) 


It is of interest to compare Equations (5.11) 
and (3.12) which represent the inplane equilibrium equations 
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with the ne of the Mushtnri-Vlanov approxi ation of thin 
shell theory specialised for the care of a shell of rtvcli- 

tion co that we let iL « r* . « r . a » <i> and c =6 

• 1 v 2 o 1 g 

(See Section 17, p.84 of Ref.fc)), This perraitc the ucc of 
the auxiliary function (stress function) 2?((},G) of the 
Mushtari-Vlaoov approximation co that the stress resultants 

^H» ^ 2 '} 11 \2 RrG de * illc d in terras of F(C|,G) ac 

follows [2] : 

M ‘ Cot * r 

N11 - -71 F ee — F > + 


1 a 1 1 


It follows from the above that 

C>v n 2; ) = - v 2 F - - — *] + F - eo } 

Since r Qt r,j and r^ are independant of tho 8 
co-ordinate in the case of shells of revolution, if f 
H 22 and N „ as defined by Equations (3© 15) arc introduced 

into tha inplane equilibrium equations (3.11) and (3.12), 
the resulting equations will be: 


(3.15) 


(3.16) 


ih;(T F ’*) = 0 
rir(r F >O +0 - 0 
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The above tv.o equations are satisfied exactly for phell3 

of revolution of zero Gaussian curvature (i.c. -- - « o » 

r r 
1 2 

as in the croc of circular cylindrical or conical shells 
provided the applied load components $ and 0 are zgjo. 
In the caoe of other shell geometries, the last two equa- 
tions will not be satisfied exactly and in such casco, the 
error involved particularly when .$ and © aro zero, will 
be assumed to be negligible. 

It is worth noting that if in Equations (3.15)r 
the function F is replaced by w, tho resulting expres- 
sions becoino that of ^ (- ^ r> ) as given by 

equations (3*5). From Equations (3.15) it is seen that 



where the operator 15 stands f( r 


(3 


DO - ) 


j_ r_a_rjo_ <K~n i a*c~o j 
r,r 0 l d, Lr,r 2 a#» J + r 0 "de 1 i 

(3 


Similarly, from Equations (3.5). it follows that 
(£" t iSll.) r -D (W) 

Because of Equations (3.15) and (3*17), Equa- 
tion (3.13) which involves the .-‘unctions w((},Q) and F(Cj,)) 
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can be witten tin: 




D ,.w - w, f ) . -I (±u a w,.)] 




3|*J_ 

a*> .v, 



(3.19) 


The last equation involving the functions W and F 
is quite complicated because of tho presence of the terms 
within the brace B on the r. ght hand side of this equa- 
tion. These terms multiplied by D(1 - y ) can be shown ;o 

stem from tho integral \ ( X . ft — K* ) dA appoar- 

J A 11 22 12 

ing in the expression for bending strain energy U 2 an given 

by Equation (2.8), If the integial mentioned above is neglec- 
ted so i-hat the bending strain energy is approximated by 

(f)J A (x n + x 2i^ dA =(T)I A (? 2 w) 2 dA , . th*i, 

a considerable simplification of Equation (3.19) results. 
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Consequently, the following simplified version of the coi- 
tion. of equilibriu.a in the ^ coordinate direction is 
obtained. 



( 3 . 20 ) 


In the last equation, ill terms on it a left hard 
side except the first two terms involve coupling between 
W end F and hence this equation is coupled# If these terms 
involvi:ig coupling are neglected, the resulting linear equa- 
tion, namely, D W D(F) = Z represents the equilJbirum 

equation for a shell of revolution according to the Uushlrri- 
Vlasov approximation of thin shell theory [2] • 


Equation (5.20) can still further be simplified by 
introducing the inplane equilibrium equations (3.11) and 
(3.12) //herein the loading terms <[> and 0 are neglect(d# 
Thus, from Equation (3.20), the following equation is 
obtained : 


DV 4 vv . 0(O-*h A(i W „)- K a [^1 W „ 

wJ = z 

r„r, L ,e * r a 0 J 


(5.21) 
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Using the ftxprefcDiono for find >C An c.c given by 

11 22 12 

Equations (3«b)# the last equaMnn can nine be written ii* 

the form 

DV*Vy-5(F) + N„ ?t,| + N 22 ?l 22 + ?.N, 2 X, 2 =■ Z 

I 

(3.22) 

It should be noted that Equation (3.22) is alee 
the governing equation for the buckling of a shell of re*\olu- 
tion provided it is recalled tliat in buckling problems, ^ 
and U 2 2 are compressive (and similarly for U ) o° that lor 

buckling problems Equations (3.22) becomes 


D(v*w ) -(—+ - N„ Jt„ -N 2 iXj,-2N, 2 = Z 


or, by means of Equations (3*5), 






(3.23) 


(3.24) 


In the above equation. IT . N and N represent 

1 1 22 12 

the stress resultants (acting on a typical shell element) 
indtieed by the applied load(3) *;liat give rise to tilling* 
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c 


Their values just prior to buo):lJng arc determined iron tie 
inplanc equilibrium equations, namely, Equations (3.11) 
and (3.12). 

Because of Equation (3 « 15), Equation (3.21) 
represents one equation involving; the two function'. XI cn< l F. 
A second equation involving XI and F is obtained from th) 
compatibility equation. If in Equations (3*5 )» the nonlinear 
terms are dropped, the resulting expressions are exactly 
those of tho Munhtari-Vlacov approximation specialised fo:? 
shells of revolution. The compatibility equation that is 
made use of in the LIushtarifVlasov approximation specialised 
for she] Is of revolution can be shown to be [2] : 


y y (l‘± , ttMV a_ril.fr c V- 1/1 Ikr '\ 
° '\y 2 r, ) a^Lr, d<t> Cr ° d <t> e ') 


Y t 3 2 €| m „ . BfO 


3 a u) 


= 0 


(3.25) 


If we linearise the expressions for 6^ , 6^ and to 
as given by Equation (3.5) and introduco the resulting 
expressions into the left hand side of Equation (3.25 )y 
the result i3 


The u end v dependant terms ariie because in writing the 


1 f 

V. 


: % i 


expressions for curvatures % , % and % r no dependant 

1 < 22 1 2 

solely on the displacement coispcnent w($, G ), [ nan Bqu.^ion 

(3.5)] , an approximation is involved and thin approximation 
is also Mfijdc! use of in tho Mushticri-Vlunov approximation. If 
tho u and v dependant terms in the curvature expressions had 
been kept, the tcrrac - ain $) + would not 

uy r 2 

arise. In what fo.V v,?s, those terms depending upon u and v 
will be neglected. In this approximate sense, the right hand 
side of Equation (5*25), upon introducing the relations lor 
e 1t and (0 as given by Equation (3.5 )» will become 

[- r„ r, c tw) + r 0 r, + ^p.)] 

depending upon W and its derivatives • 


plus nonlinear terms 


From Equations (3.9) it folloY/s that 


(n„->’n 22 ) 


•I > 


£ 2 = ^(N»-r^)j « = -^- 


(3.26) 


Equations (5.15) and (5.26) yield 


c , , I c Coti y 3 / 1 c N 

Eh£| r* F, w r,r 2 F >* + r, s<j>\Y, E, v 


Eh6 -"T] 


ycot d c 
91 r, r 2 ^ 




EhB- 20.«P T F, M -^F„] 
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Introducing Equation (5*27) into the left hard 
side of Equation (5.25) tho following cxnreosion is obt c?ncd: 


Eh Id^Lr, 3(5 L r, a^v»i '*/ r 0 r, ’M J 


3 pCos^ c Cos 2 $ r ^CoS(5 3/1 V 

+ 3(5L’ _ tf - + IffT •* ~ ~T ' _ W ’+)■ 


1L' + S±t( F ',-Z L(Lf V 

r 0 D? »•**• /,*•> W w ./ r, d*Vr, + c6 /_ 




-20+3')^ (j". 4-ee) " ( r< ?^) F>e »] } ^ ^ C‘ 


^ %YL ' 


The above expression can be shovn to be equal to 


. Y Y (Xtl . #** 

r ° ‘ l"T + “^" 
\ >2 ' 


where , 


y4 p = V 2 ( V 2 F ) and the liaplacian operator v' 


is defined as 


r a pb ac -oi 

Id^Lq d<f> J \r 0 y ae 2 J. 


Introducing 6 1f 0 2 and cO as given by Equation 
(3.5) I’.nto the left hand side oT Equation (3.25) and 
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neglect ii.g, an mentioned onrllcr, the termc 


[-.VS- *">)-« 

expression i 0 obtained: 


» the foiiov/inr 


-ror, SM ♦ {(&)[£ 

f 

+ {($(£ w - ») tA' 7. w fa <i v fa T f (% «. a) ~(f) 
+ ^ M ' v ’ * ) " K •») fa *) ft (fi . ;) 

(Coj (!) [(— W, f ) — (y \N,^j +-L (\V,t)lW, $e ) 

+ -^ (W, e ) CW, dji)] ]• + Ke,r, 


On equating expressions (3.28) and (3.29) tho 
follov/in,; compatibility equation ia obtained: 

Et (y4f) + 5<vv,= ^{(t)LVi w -')] 2 
+ U fa- «) fa’ fa fa fa fa, fa ] } f 


0.30) 


In the above equation, the right hand aide ccnoiGtn 

I 

of nonlinear terms and if they are neglected, the resulting 
equation is that corresponding to the Hushtari-Vlasov approxi- 
mation of thin shell theory [2] .. 

Equation (3*20) is too complicated and a simplifi dc.- 
tion of it is essential if any solution in to be found for on 
actual problem. In ouch a simplification - to be justified 
later - all the nonlinear terms ccntaincd within the brace 
♦on the right hand side of Equation (3*30) are neglected, Cf 
in the remaining two nonlinear terms, we further approximate 
and write Yt T ,G) ~ ■“£“"* Y/,Gj) which means neglecting 

W,0 ,dr n 1 1 

the term — 7 ? in comparison to (jrj \7,Q<p) f 

Equation (3»30) will finally simplify to: 


V 4 F +E 


h [ 5 H “ 


(3.31 ) 


The basis of the simplification mentioned above 
requires justification. To this end, the compatibility 
equation (3*23) in re-written as: 

_LLLjL Ci e 0 jL{±l dr °) i. 

a<#.Lr, ° 2 \ a^Vr, dyi/ v r c ao* 

_ L (QsS 3(3 'fj’L = _ r r (Jk. , 2M 

ib Vd i) ao d4>de 0 1 V.’r 2 r, ) 

J . 

If in the above equation, terms involving products with 
are neglected, the resulting simplified compatibility 
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equation is: 


r 0 c n $ i i 3*^ 9 a to 
r, dj6 l r D ae a z? do 


/ ?C|| 

« -r 0 r, (— + ---) 
\ I2 r, / 


If in tie above equation, the nonlinear part of the expres- 
sions fcr e 1t e 2 and co as given by Equation ( 3 . 5 ) arc 

introduced, the left hand side of the last equation yields: 

On expanding the above expression and neglecting terms involv~ 
ing products with ^ or » we will finally obtain 


the expression, 


p Q r, [( W,e< ^) WW] 


T?V\ 

The abovs expression multiplied by is the right hand 

o 1 1 

side of Equation (3.31 ). 

Finally then, we have Equations (3.21 ) and (3.31 ) 
as the governing equations (for the static case) for the 
nonlinear analysis of thin elastic shells of revolution. 
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They are a pair of coupled, nonlinear partial differential 
equations* If the nonlinear terra are neglected, Equations 
(3.21) and (3*51) reduce to tho governing equations of 
Mushtari-Vlasov approximation of thin ehell theory specir.** 
Used for the case of shells of revolution. 

It is worth mentioning Ifcat Equations (3.21) and 
(3*51) also represent tho governing equations of the problem 
of poctbucklir.g behavior of thin elastic ehelle of revolution 
provided the signs in front of tho terms involving the stress 
resultants and N defined in terras of the strcsB 

function F($>,0) by Equation (3. 15) are changed from minus 
to plus. Thus Equation (3.24) holds good; only, the stress 
resultants are to be determined from Equations (3.31 ) ar.d 

(3.15). 

For the analyses of nonlinear flexural vibration 
problenf of thin elastic shells of revolution, Equation (3*31 ) 
remains unchanged as it represents the compatibility equation. 
However, in Equation (3.21 ) Z has to be replaced by 

r 2, 

e, t) - p ] where, p is tho mass per unit 

area of the shell and t is time. For free vibration 

problems, Z ( 0. Q, t) ■ 0. 

. applied 

It is worthwhile summarising the governing equations 
for free undamped nonlinear flcx'iral vibration of thin elastic 
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she 11 o of revolution* There equations (in orthogonal 
co-ordinates <J> and 0 ) involving the normal displacement 
function '({>,0,t) and the auxiliary funotion (stress fun:s- 

tion) G.t) are:- 


V*F - Eh {-6cw) + (^w,J 2 -( T y w ,„)} 


(3.32) 


(compatibility equation) 


pw t D(^W)»(5(F)t^fiw,^ 

+ N « IrSr K?i + ¥ w > J + 2Ni2 [« w ’°* ■ ‘T i & w ’®)l J 


(3.33) 


(Equilibrium equation) 


The 8tress resultants N« < , and II are defined 

11 22 12 

in terms of the stress function F by equation (5.15)* 
However, in view of the approximations involved in Equations 
(3.32) and (3.33)» these stress resultants can be approxi- 
mately expressed by 

N|| ” ‘jl ^ 00 > ^22 * _ *j^2 M ') N I2 25 jT <f>f 

Also, the two terms in Equation (3.33) involving 
Cos <f> 


os l /dfcN 

r 0 “ r 0 r, Vd^/ 


can he dropped. 


HI 
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The biharironic operator is defined by v r = v^- 

J> 

where, the Luplacian operator v‘ in defined by Equation 
(3«20a)* But in view of the approximation employed in obtain- 
ing Equations (3*32), (3.33)» .and oo on, can be 

approximated as 


V 2 (••• ) 


n aHO _j_ tfc-n 

U 2 a? + t i 2 ae 2 J 


The above equations can be readily specialised for the 
caoe of any shell of revolution; for example, circular cylin- 
drical, conical, ogival, torus, spherical shell and eo on. 


With do^ =» r d{) and do^ « r 0 dO, Equations (3*32) 
and (3.33) can also be written at: 

* - [JL-'JL) C 

(3.34 ) 


’‘'■“t-'-yar-oe©} 


and, 



t !)(V 4 W) = -|d (F) + Nii 


3*W ri (d r o\ 3w 3 S W‘ 

3s, 1 ri l r 0 Ids,/ 5s, + as£. 



a a \v 

05, d >2 



(3.33) 


where, 




and frera Equation (3 .28a) 



(3.36) 
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since terms involving products of 


8r 

0C. 


o 


ecu be neglected 5.n 

dr. 


1 ui o 

which case the two terms proportional to (dr* Z7~~) in 

'o Uw i 

Equation (3*55) can also be neglected* 


It is of interest to compare Equations (3*34) aid 

(3*35) with Equations (2*34) and (2*35) applicable for the 

case of thin flat platen. If in Equation (3*34), the terr 

D(W) in dropped, the resulting equation, formally at any 

* rate, coincides with Equation (2*24). Similarly, in the 

absence of the term D(F) and the two terms involving the 

products of Equation (3*35) reduces to Equation 

r o QB 1 

(2*35). Thus the influence of the shell geometry is brou£ht 
in essentially by the terms D(W ) and 13(F), where the 
operator 1) is defined by Equation (3*10)* 


As outlined at the end of the last section dealing 
with thii flat plates, here also two methode of solution of 
Equations (3.32) end (3*33) are possible. One is to aesune 
a deflection function Yf that satisfies all the boundary 
conditions of tho problem and based on the assumed YV, find 
the stress f miction F that satisfies the compatibility 
equation (3*32). In general, the assumed \! and the corres- 
ponding F will not satisfy exactly the equilibrium equation 
(3*33). The Galer kin's method in employed to obtain an 
approximate solution [ 3 J » [4] • 
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The alternate Method in br<sctl cn writing the 
inplane equilibrium Equations (3 .11) and (3.12) in terms of 
the displacement components u< v and \i oo that from the 
resulting two equations, we can try to find u end v for p.n 
assumed V7. For example, if in the inplene equilibrium 
equation (3.11), Equations (3.S) tire introduced and the median 
surface strains 6^, and defined in terms of u,v and w 

by Equations (3*5) are wade use of, the resulting equation 
can be shown to be: 
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Similarly, stnrting out witli Equation (5.12), another 
equation in terr43 of u, v and w can be written. The 
right hand sides of these two equations involve VT and 
its derivatives so that from th;ce two equations, we can 
try to find u and v for an assumed YJ. This method ha 3 
been used in the solution of the problem of nonlinear 
flexural vibration of oval cylindrical shells [ 5 ] * 
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Having derived the governing equations for the 1 ii'gc 
amplitude (nonlinear ) flexural vibration of thin elastic flat 
plates using curvilinear orthogonal co-ordinates and 
as well as of thin clastic shell j of revolution, it is possible 
to generalise and write down the governing equations for the 
large amplitude flexural vibration of anv thin e?..astic shell. 
From these very general equation!, the governing equations 
applicatle for shells of revolution or thin flat plates can be 
obtained as special cases. 

For any thin elastic shell, let a typical point on the 
median surface be P. At ?, the directions of the principal 
radii of curvatures rhich are orthogonal and the orthogonal 
co-ordinate directions and a, ? on the shell median surface 

are assumed to coincide. The elamental length ds on the 
shell ra ;dian surface ia given by 

(ds) 2 - (A,doCf) 2 + (A 2 doc z ) z ' 

= (ds,>*+(cfs 2 )* (4- 1 ) 

where, and A^ arc the lame' parameters and arc in general 
functions of a-j and a 2 * The principal radii of curvatures 
at a typical point ouch as P are denoted by r,j and i % 2 and 

they too are in general functions of and a 2 . The two 
Lame' piraraeters Aj and A 2 and the two principal radii of 


PADS® 

ftUAun 
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curvature 3 at a typical point ouch as ? satisfy the three 
conditions luiown aa the Gauso-Codai;?.! conditions, nane'jy , 


a f^- \ _ J_ fslhi. jL_ f A| V j_ 

doc, \ \'i) T, wet,/ 9a 2 \ r, ) ti \d«2/ 

(Conditions of Codnnzi) 


and, 


"a ^ i aAA + 3_ /j_ 
a«., sA| 3ot,/ da 2 VA2 



(Condition of Gauss) 


(4.2) 


(4.3) 


Tho strain at a general point (a^, a r/ , y ) can he 
Y/ritten in terns of the six median surface deforuiationnl qi anti- 
ties 6.J, e^, to t #i<j» ^2 311(1 ^12 33 Equation (3.4 )» 

where [l ] , 

_ J_ 3u. _T5_ /3A,\ W J_M _ 3lV*1 2 

^ A, 3a, A|A£ \3 oc2 ' H 2l A, 3aji 


I dv t U f dA 2.) { _VV JM 3W1 2 

A 2 A, A 2 \^0L^J Y>i 2 L A 2 

fi u. /^3A,\ au- _ tf_ f aA a \ n /awvaiv vn 

[A, da., A, A 2 '3c*2' A 2 3«2 A,A 2 \da,y [A,A 2 \3oL|A3oC2/Ji 



(4.4) 
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The last eet of equations specialise to thane 
applicable for the case of thin flat plates , namely, Equ - 
tions (2.5) and (2.4), provided It is recalled that the rr.dii 
of curvatures (r^ and r^ ) of the flat plate median surface are 
infinity . With •* «• and r 2 ■* «, Equations (4.4) reduce 
to Equations (2.5) and (2.4). Furthermore , with A 1 « r 1 and 

A 2 a r o as wel1 as a l * $ ^ a 2 c Equations (4.4 ) 

specialise to those applicable iur the case of thin shells of 
revolution, namely, Equations (5.5). It should also be 
remembered that for shells of re/olution, r^ and r Q (and 
henco r 2 «* are functions of $ only and not of 0 

coord incite* 

From Equations (4.4), it follows that 


(±0 

?JV l 

+±.p 

^ aw ~[| 


Lacd, <■ 

y 3a, j 

aa2.Lv/ 

i 2 / aaj J 

(4.5) 


Also, 




, J_ /in/M dvr 1 i_fi/ r M§w‘n 

A| A 2 . L ^oj.|L A| / Bcx»| J 8a2L^| v A 2 / d0C2*J J (4.6) 


It is to be noted that in the case of a spherical shell, 
r^ = r 2 35 R Y/here R is the rriius of the spherical shell 
bo that the operator D (•••)" ^ ^ 2 (- • • ) 

Using Equations (4.4), the energy method and ntepo 
analogous to those outlined in tne last section leading to 
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Equations (3.32) and (;>*33)» the tfoverning equations for the 
free, undamped, lar^o amplitude (nonlinear) flexural vibra- 
tion of thin clastic rhcllc con he readily obtained in 1 iritis 
of tho normal displacement function V(c < ,, t) a' d tic 

auxiliary function (otreay function) F(tt^ , cr 0 , t). These 
equations can be cho\m to be:- 


( 


r % ( I a 2 w \ 2 f \ a 2 wv I a a W\) 


(Compatibility equation) 


^ Nil 8 / I 8W \ ^ Nji 8 / 1 8W 


At d0C|\Aj dtt|/ d0ti'‘A2 d((2 


AjA;. dotfd&2 


n aw\ 

VA2 8Kj/ 

5} 


( Equilibr ium e quat ion ) 


y/here , 



r 1 a r i bf ^ i /8 a, y i aF \i 

lA| 3 a,\A| Dcl 1 J A,A2 \8a2AA2 dca^JJ 


and, 

i__ r d*F f aA ' Y' 3 H W 1 aF T 1 
l2 " + A,Aa LdA,0«2 \.3«. 2 Aa, 3a.,/ I 3oL, A A 2 da 2 /J 


(4.9) 
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However, in *rJLew of the approximations involved in Equ^V.ons 
(4.7) and (4.0) which arc tho r.?.ir.e no thoso involved in 
obtaining Equations (5.52) and (3.33 )» the stress result .vrtr: 

ORIGINAL PAGB 1H 
04-' POOR QUALfCff 

Kj . ~ ' : . (4.9a) 

* 2- A| Aj. dOC.|dOl2 

The approximation involves neglect of thoso terras involving 
products of the gradients of the; Lame* parameters, namely 



It follows from Equations (4.9) that 

^ N 11 + N 22^ ** “ V 2 F 011(1 **■ ( W 22/ r 2^J “ " 

A A p O 

The biharmonic operator V f is defined by V 3 V* . V 

p 

where, the Laplacian operator v ic gi ven by Equation (4.5). 

Equations (4.7) and (4.8) can be readily specialised 
for tho case of thin flat plates by letting « and 

r^ - 00 oo that L(Yi r ) • 0 and 15(E) * 0 and consequently, 
Equations (4.7) and (4.9) reduce to Equations (2.50) and (2.51 ). 
These equations also specialise to Equations (5*52) and (5.35) 
applicable to thin elastic shells of revolution by letting 

A 2 3 r o* a i ” ^ ancl a 2 = ^ should be recalled 


can bo approximated by: 

. _ i a_ n ar\ . 

11 ^ a 2 d<*i \Aj> act 2 / 9 
K, I B f\ aF \ 

21 ~ a, aoc, V a, aoc, / 9 

A A J 
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that wilhin the frame work of the approximations involved 

in obtaining Equations (3.32) and (3.53), the two terms 

proportional to £ 2 ?- .^ « ----- fo^*) should be droijped, as 

r o 1 o 1 1 

mentioned in the last section. 

A more accurate equation of equilibrium than Equation 
(4.3) cfJi be shown to be: 

pV + DCv^W) = [d (F) - Nl„ - N 22 K 2 z'2H ll ?C i ^\ 

1 “ J (4.10) 

where, and are given by Equation (4.4). 

For the problem of buckling of the shell (either the 
linear buckling problem or the postbuckling problem), Equa- 
tion (4.10) without the inertia term pw is applicable p rovid ed 
the signs in front of and are changed from minus 

to plus so that for shell buckling, the governing equation 
becomes. 

D (V 4 W) + D (F)-N||?t,| — N 22 - 2 N| 2 %i 2 - Z 

(4.11 ) 

where, D(F) « - [(N^/r^) + ( W 22^ r 2^1 *nd z represents 

the component of the applied surface load in the direction 
of the coordinate. For the linear buckling problem, ;he 
values c-f the stress resultants lfj ^ and just pr .or 

to buckling are determined from the linearised version of the 
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coinpat ibility equation (4*7) wherein the term l)(YO ic 
dropped. The relationship between the stress Tunc t ion F and 
the stress resultants are define i by Equation (4.9 )• I or the 
postbuchling problem, the stress resultants appearing in Eq na- 
tion (4.11) are determined from the compatibility equation 
(4.7) which is nonlinear. 

f 

With ds 1 « A.j da,j and dc., “ da^, Equations (4.7) 
and (4.N) can be written as: 




(4.12) 


and, 


{ O, 

DCF)+H„ + N 2 i 
ds i 


In the above equations, v « 
equation (4.5), 


d\-J 

ds/ 

2 


+ 2 N|2 


3 2 

35, d5£ 


1 


and from 


V 2 0 •) 


J_L 

, 3C_)1 

1^>2. 55| L 

' 2 ds, J 

p a c-> 

a a c— >i 

~ Lasf- 

3sJ J 


I 


A, 3^2 




(4.13) 


(4.14) 


since l ent s involving products >£ the (prr.dientB of the Lime 7 

6A„ 3A, 


parameters, namely, and 

OS.j 


y/u * 

•r— are neglected in the 
0s 2 


approximation employed here. 


With N 




d r f 


ll - 


a si 


N 


22 


' tr? a " d 


HTF 


as, as* 


Equation (4.13) can be written a3 


p»*D<7«W) = T DCF) ^ 

1 1 \dsU asf \ds,V a$| \a5,a5 2 / ac,3s 2 J 


We now define an operator (following Vol'rair) 

OCW,f ) = F 5l5i+ (W^) F 5 | s, -2fs,s 2 Ws, 5z 


bo that we get, 

pWtD CV^W)^ |d(F) - 0 CW,F)j 


/ 


which i£ a neat, chort way of writing Equation (4.13). 
Since, 


OCw,w) = 2 [(.4 , s ,Ho 2 " fan) 1 ] 


we can vrite Equation (4.12) as 

V^F = Eh|-Dtw) -i 0 CW,W)j 

Equal iono (4.12) and (4.13) hold good for any thin 
elastic shell and it is readily verified that from these -qua- 

0 

tions, Equations (2.34) and (2.35) applicable for thin flat 
plates and Equations (3*34) .and (3.35) applicable for thin 
elastic shells of revolution can bo obtained as special crises. 
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Section 5« 




In this section, on tb* basis of an assumed vibration 
mode which is taken as the product of a function of time e.nd 
function of the space variables, the modal equation is deter- 
mined for the £eneral case of thin elaotic shells discuss* & ir* 
the last section. This modal equation can then be specialised 
for shells of revolution or flat plates or beams. The ch*»rac~ 
teristic features of the modal equation as applicable to rhin 
* shells cn the one hand and thin flat plates and straight beams 
on the ether hand are discussed. 

Let, 


W(s,,s 2 ,t) = B (t) W (s, , s 2 ) 


(5.1) 


where E(t) is a generalised co-ordinate and , s„ ) is 

the mode shape assumed to satisfy all the boundary condit-ons 
(both geometric and natural). Introducing Equation (5.1) in 
the compatibility equation (4.12), v/e get 


V 4 F = Eli|-B«) &(\V) + e? 


17 \ 

( a 9 ' w\ ( a a wl 

Lvas, Dsz / 

WAasi/l 


(5.2) 


where 



D-OT>={- 


a aw i + j_ a r/ A t \ awl] 

5sj a y<2_) 36j- a i as 2 L xj dsz. J 


(5.3) 


Equation (5.2) is satisfied provided F is assumed as 
F Os , Sj , t ) = B (t) F (S.jSr) + B 2 P (s,,sz) 


(5.4) 
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So that, 


V 4 F =• -Eh[D(\V)] 


and 


V'*F 


Ehff — f 

Lvas, 0 s 2 / \d$j* 


s )@)] 


(5.5) 


(5.6) 


If the problem were lir.eor so that Equation (4.12) 

B 

ifl linear, then, 3? « 0; that it, FCs^Cpjt) = E(t) F(s.,c 0 ). 
The equilibrium equaticn (4.13) can be written as 


L(Y/,F) « 0 


(5.7) 


where , 


L(W,I : ) = | pW +1)^^) - D(F)— Kl, 


„ -N„ 2N, 2 

as* 22 as* 2 


b 2 vj 


gNv \ 

Ac . as* 1 




Uoing the Galcr kin's method, 

4 [l(W, F)]Sw. ds = 0 

where, dS a A^Ag da^ « dc, ds^ and 67/ a (6B) W 
bo that Equation (5«9) becocer>: 


J [l(w,f)J W ds = 0 

S 


(5.10) 


By meane of Equations (5.1). (5-4) and (4.Su), Equation (5.C) 
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yields: 


L(W,F)=j;pS)B+ b[dv‘c®-&cf>] *“*[©(S)*(SXS') 


-€k) tsj- 6 <" ] > 43)(0MS)@ 


_;.f£L4 (*'* \] 1 

'dSj dZ>2. '9Sj dS% J . j 


(5.11) 


It is helpful to define the follov/ing integrals: 


I, = P I (W ) 2 ds. 


L Z = vf s C'^ W) wds. 


Ij= f [DCFj] w ds 


WM 



9 ** 


d z F 

d^W 

a 1 f 

a*F 

a^'vv "] 

ds| 

dsf 

as* csf as,3s2 

ds, as* 

a^F 

/w 

-2^ - 

afo ■) _ 

Vw 

ds 

ds* 

as| 

dS,d:i2 3S,352J 


W ds 


(5.12) 




r 

\ 
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In terms of the above quantities, Equation (5.10) can be 
written as: 

j! 

ii 

I, 8 +(lj — Ij) B -t- = 0 (5.13) 

The solution to the above nonlinear ordinary differen- 
tial equation involving B(t) helps determine tho amplitude 
versus frequency relation in tho case of thin clastic shells. 

If the shell v/ere undergoing bs ill amplitude flexural vibra- 
tion (that is, the linear probljn), the associated frequency 
equaticn can bo determined from the linearised version of tho 
modal Iquation (5.13) bo that the governing linear differen- 
tial equation becomes 


1,6 +U 2 .-I 3 ) B = 0 


(5.14) 


Assuming solution of the form 3i(t) a B Q coo cot, 
quency go as determined from Equation (5.14) is 

i 


60 


tVr 


the fre- 


(5.15) 


Tho nodal equation (5.13) can be readily speclaliccd 
for the case of shells of revolution for which A 1 K r -j » 

A 2 ** r o» a 1 3 $ an< * tt 2 8 6 so that do^ * A^ dc^ ■ r 1 d$ 


and de 5 = Ag da^ » r^ dG. Also, for a shell of revolution, 
the elemental area do a r Q r 1 d$> dG. 


- 
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It in of interest to obtain the modal equation for 
the (nice of thin elastic flat platen. This can be obtained 
in a manner analogous to that used to obtain Equation (5.13), 
and staining out with Equations (2.34) and (2.35). The 
required equation can also be obtained by specialising Equa- 
tion (5.13) for the case of flat plates for which the radii 
of curvature of the median surface r^ and r^ are infinity. 

When r^ - «> and *• », it is obvious from Equation (5.3) 
that , 

D(w) = 0 and D (F) =0 (5.16) 

Thus the compatibility Equation (2.34) is satisfied by 

F (s.,Si,t)= B z f(s,,s 2 ) / (5.17) 

comparing the above with Equation (5*4), it is seen that 
for the case of flat plates, 

F C s, , s 2 ) * 0 (5.18) 

Because of Equations (5.16) and (5.18), it is clear from 
Equations (5.12) that for the case of flat plates, 

1 3 s 0 ; I 4 * 0 (5.19) 

Consequmtly Equation (5.13) simplifies to 


I, B + I z B + I 5 B 5 = 0 


( 5 . 20 ) 
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The above modal equation which is a nonlinear ox di- 
nary -h inferential equation helps determine the anplitude 
versus frequency relationship for a thin clastic flat picte 
undergoing large anplitude flexural vibration. On conpajing 
the last equation with Equation (5.13), it is seen that in the 

case of flat plates, the nonlinear term proportional to B 2 is 

# • 

absent. If the last equation is linearised, the governirg 
equation for the small anplitude flexural vibration of a thin 
elastic flat plate is obtained, from which, the frequency GO 
(rad/second) can be determined as 


» =An, 


(5.21) 


The nodal equation in the case of straight beams 
undergoing large amplitude flexural vibration can also bo 
shown to be given by Equation (5.20). For this, it is 
easier to work with the governing equation for the large 
amplitude flexural vibration of straight beams, namely, 


.. d*w 

P w + E1 l^r - N <> 


3x* 


= o 


(5.22) 


where, « Ea and V is the beam cross sectional urea. 

From the equilibrium equation in the x- coordinate 

dG. 

direction, it can be shown that - — L *0 so that G *= C » Con- 

dx lo 

stant. Thus, « Ea 0 o v/hero G^ £u, x + l(w, x ) 2 J. It can be 


t 


( 

V. 



shown that 


<) 

If Cu)' = 0 , 


c«= ~S. wf K ) clx 

=z H +k£(*M- 

c o*r/X( w >*f £Us » that N « = ff IWt/ ds ' 


Assuming W(x,t) « A(t). \7(x) and using the Galerkin'c method 
as before, the modal equation can bo ehovm to be 


J, A + j£ A + J 3 A 5 = 0 


(5.23) 



(5.24) 


It is of interest to note that the modal equation 
[see Equation (5.13) or ( 5 . 20)1 can also be obtained by apply- 
ing the Lagrange's equation of motion, namely, 


dt 



0 


(5.25) 


where, T is the kinetic energy of the system, U represents 
the strain energy, B is the generalised coordinate and there- 
fore, B is the generalised velocity. With, 

t = 1j C<-^ ds 

where W is given by Equation (5.1 )» we get 
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In the above, I s is given by Equation (5.12) so tliat the 
above represents the first tern of Equation (5.13)* Similarly 
^ the other tenia of Equation (5.^3) or (5.20) can be sho.vr. to 
stein from the second tern of Lagrange’s equation given b} 
Equation (5.25). 

Thus, the modal equation from which the frequency 
versus amplitude relationship 1* determined in the case of the 
large amplitude flexural vibration of straight beams and 
« plates is of the same type [sec Equation (5.20)] whereas, in 
the cao 3 of thin elastic shells the governing equation [nee 
Equation (5*13)] has an additional term proportional to B . 

It is worth investigating whether this circumstance has any- 
thing to do with the fact that in the existing literature, 
most of the available solutions for the case of thin shells 
are found to give rise to nonlinearity of the Softening 
type" (that is, frequency decreases with increasing ampli- 
tude) vhereas in the case of flat platen of different geo- 
metries the nonlinearity is alw<iy3 found to be of the 'harden- 
ing type" (that is, frequency increases with increasing 
amplitude). In fact, based on the survey paper [fc] and [v] , 
it io i-oosible to draw the following general conclusions 
regard! ng the type of nonlinearity exhibited by the different 
elastic bodies. 
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(A) Straight beams - . Hardening typo 

( 1 ) Rectangular 

(2 ) Circular 

(3) Elliptical 

(4) Triangular 

(5) Skew 

(6) Rectangular 
Sandv/ich type 

(C .) Rings - - — - Softening typ •* 


(B) Plates 


> 


-- Hardening typo 


(D > Shells! 


(1) Circular cylindrical 

(2) Oval cylindrical Softening type 

(3) Spherical 


The above oumory of results raises the question as 
to why in the case of straight beams and flat plates (irrespec- 
tive of plate geometry), the nonlinearity is always of the 
hardening type whereas in the case of bodies like rings and 
shells the undeformed median surface of which have curvature 
to begin v/ith, the nonlinearity is of the ooftening type. 


Equations (5.13) and (5.20) can be written formally as 


d 2 0 

dt 4 


+ ale + be 2 -i ce 3 = o 


where, a Q ,b and c are real numbers and 0 « 0(t). 


(5.26) 

It will 


* In the case of cylindrical shells, it is assumed that 
the ratio of the axial to circumferential half-waves 
of t le vibration node ir> small. If thi3 ratio i3 large, 
the nonlinearity is found to bo of the hardening type. 
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be easily recognised that the abrvo equation also represents 
the dynamic equation of a mass-spring combination whore the 
b spring restoring force is nonlinear and is defined by a cubic 
expression of the form (a Q 0 + bO + c0 p ). 


,d0 


Multiply each toria of Equation (5.26) by dt ) and 

integrate to obtain 


(0) 2 =C o -[a« 0 a +|-be 3 +|*0 4 


where C 0 a constant of integration to be determined from 
initial conditions on 0 and 0. From the last equation, 
we can .finally obtain 

.0 


t- r£= 

e # Soce) 


(5.27) 


where C> (0) = [c o - (aj e 2 + | b O 5 + | € 4 ) 


is a quoriic 


in 0 and G q w 0(t = 0). The integral appearing in Equa- 
tion (5.27) is an Elliptic integral of the first kind. The 
nature of 0 versus t curve will depend upon the constants 


2 

a Q , b, c and d 0 « For the flat j.lato problem the governing 
equation of which is Equation (5.20), the coefficient b « 0 
in Equation (5.26) so that Q(0) appearing in Equation (5«27) 
becomes a biquadratic. This is true aleo of the problem 
involving straight booms [see Equation (5.25)] . 


When b =* 0 in Equation (5.26), the solution to the 
resulting differential equation where c > o can be chovm 
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to be (nee page 26 , Ref.Q) 


t = 


(a* t Co?)'* F0 ' ,9 ' ) 


where, 


C r 6^ 

F1 ^> -i/rSE 




(5.2B) 


(5.29) 


Je the elliptic Integral of the firat kind and 


K- 


where a. = 


C© c 


20-ta) ~ a* ( 5 . 30 ) 

obviously X 2 < 1 . The initial conditions are assumed to be 
G(t *= 0) =* e 0 and 0(t = 0) « 0. It can also be sho\rn [g] 
that the period of oscillation of the dynamical system for 
which the differential equation (5.26) with b « 0 applies 
is given by 

X - 4P 

•" a 6 yn^r (5.31) 

where, l‘(A ,~) is the complete oLliptic integral of the first 

ki).d \vi ;h modulus X • Furtherw ire , it is known that the 

o 

period Z decreases with increase in a Q , <J, 0 Q either indi- 
vidually or collectively [e] • fhen G =* 0, we have the 
linear problem so that 


( t ) -■ 7-* - 2TT 

iwifow- T SLo 


Thun, 


(£W 


2 


'LL 


yr?a (5.32) 

The above ratio which i« equal to unity for C « 0 
can be shown to be less than unity for all c > o [o] • 
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Consequently, the associated nor -dimensional frequency 

Coo/co*) ~ Ct7t) where co io the frequency for the nonlinear 
^ problem and of is that for the linear problem, increases. 

Thus the nonlinearity is of the hardening type. For any fixed 
set of values of a 0 and c, since the period T decreases with 
increasing amplitude G o , the nondiaensional frequency (<0/io*) 

will increase with increase in the amplitude 9 Q which is 
exactly the result mentioned cat licr for the case of flat 
plates. Thus we arrive at the important conclusion that in 
the case of large amplitude flexural vibration of thin flat 
plates find straight beams for which the modal equation has been 
shown to be given by Equation (5.26) wherein b = o and 
provided c > o, the frequency increases with increasing ampli- 
tude (herdening type of nonlinearity). Note that this conclu- 
sion is of a very general nature and is not restricted by any 
particular plate geometry. 

A conclusion of a similar broad nature in the case of 
thin ehells would be ioosu interesting, if such can be found. 

To this end, we examine the modal equation (5.13) applicable 
for any thin shell where the assumed mode shape in given by 
Equation (5.1). Equation (5.13) is of the general form of 
Equation (5.26) and by comparing these two equations, we ace 
that 

a-o 55 C I 2 -U> (5.33) 

where I 2 and are both positive and it is assumed Ig> 
which con be physically justified from the frequency for the 
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linear problem as Riven by Equation (5*15). The coefficient 

1^ which dependo on the quantity E(1‘) vanishes for the cace 

)of flat plate b as cr.n be seen .from Equation (5.3). Thus the 

2 

value of the coefficient n Q for a thin elastic shell is Iocs 

than that for the case of the corresponding flat plate obtained 

by putting the principal radii o 1 curvaturo and r,, each 

2 ■ 

tending to infinity. If a Q decreases, Equation (5.26) with 
b *= o will give rice to a large value for period X which means 
smaller value for the frequency. Thus the trend of the non- 
* linearity is towards the noftenirg type which agrees with the 
result of almost all the investigations so far on thin elastic 
shells of various geometries. But thie result based on Equa- 
tion (5*26) v.'ith b » o i3 not quite conclusive because in the 

case of thin shells, the modal equation given by Equation (5.26) 

2 

has the term bG in it. With \> f o, a general solution to 
Equation (5»26) seems difficult. However, let U3 examine the 
equation 


d 2 0 2 , 2 

+a„ e + be 2 = 0 
d b 


( 5 .? 4 ) 


p 

where, a Q and b are both posilive numbers and the initial 

t 

conditions arc G(t » 0) « 0 Q and 9(t * 0) = 0. Assuming a 
solution of the form 

0 = A 0 + A, Cos COt + A 3 Cos Scot 

it can be shown (sec page 51 of Kef. 8) that 

(8H-, ICS*] 

where to is the frequency of thi nonlinear problem and 


(5.35 ) 
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CO ts is the frequency of the linear problem, iron equa- 
tion (5.35) it is readily seen that the nondir.cnsional frequency 
decreases with increasing amplitude G 0 which again imp? .Us 

softening type of nonlinearity. Thus in the case of thin shells, 

2 

the decrease in the value of the coefficient due to 1, 

_ 0 2 5 
[see Equation (5*33)] and the presence of the terra bO ir 

equation (5.26) each contributes individually to the softening 

type of nonlinearity. 

So far, we have worked with an assumed one term node 
shape given by Equation (5.1) and determined the modal equation 
for the generalised co-ordinate B(t). The nonlinearity of the 
problem makes the analysis based on an assumed mode shape con- 
sisting of a finite number of terms very difficult. Quito 

/ 

naturally, the analysis becomes more complicated as the number 
of terms in the assumed mode shape increases. All the sane, 
it is of Interest to obtain and examine the nature of the modal 
equations that result when a series is assumed for the mole 
shape. Therefore, let 

W(s,,s z ,t> = £ B,(t) w.Co.,5!) (5.30) 

where, B, ' s are the generalised coordinates and V/^3 which 
are functions of the space variables s^ and s^ are assumed 
to satisfy all the boundary conditions (both geometric ani 
natural). On introducing Equation (5*36) into the compati- 
bility aquation (4*12) and examining the resulting equation, 


i 
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it will be seen that the stress function F has to be 
of the forr.i, 

FCSi>s».t) “{eJb* ?,(»„*»> + (b/ fJ| + ET; B; Bj F* 1 

1 L J \ t j J (5.37) 

where, the three sets of functions {F^} , } ahu {F* ‘ 

v/hich are functions of only s 1 and are determined from 
the following squat ions 


. V 4 (Fi) = -Eh DCw*) 


(5.38) 


V<f ( ) . «,[<» /-(«, )C® 1 )] 


SvJ = 


(5.39) 


■ : 0 . 2 Eh [<%,>(%, HV^)] ».«> 

Because of Equation (5.36), 


(5.41 ) 


where, p takes on (one at a time) each one of the values 
that i can assume. Thus, Equation (5*9) yields 


/ [l(w,F)J W(j ds = 0 


(5.42) 


where, L(\7,F) i 3 given by Equation (5.8). If V l as given 
by Equation (5*36) and F a3 given by Equation (5.37) arc 
introduced into Equation (5.8) and recalling the definition 
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of Strega resultants in terms of ? given by liquation (/ ,9a) 
Equation (5*42) can be used to yield the following modal 
equation: 

y 

{( T ^ t A) + [P*V " Bi +(i 4 )*p^ t 
+ ih\ P (b t ) 3 + 2j [a 6 ) up (b i bj> + a 7 )ij P cef bj)] 

+ Zj >Z k ci 6 ) ijk|> ( Bi Bj B k )j = o (5.43) 

In the above, p take3 on (one at a time) each one of the 

values that i takes and J * It £ i. The coefficients -’I. ) , 

1 iP 

(l 2 )jLp p tc, are surface integrals, the integrands of which 

— - 3 . * * 

involve functions of V^*a, F/s, F^ o and p »£he exj)licit 

nature of the coefficients are rot given here, as they are 
complicated expressions. For each value thac p can take, 
we can write an equation using Fquation (9.45) so that there 
are as many equations a3 there ire p's (that is, i's). Vhio 
set of nonlinear equations are coupled in the generalised 

i 

coordinates B^s. If a one tern solution is c msxdered 
so that i «* 1 (and since j £ k £ i eo that the quantifies 
with j and k a3 subscripts do not exist), Equation (5*45) 
will reduce to Equation (5.15)* From Equation (5.45 ) iv is 
readily seen that the nonlinear terms are of two types, 
namely, of the quadratic type [hf and I3 i B.. terms] and of 
the cubic typo (B? f B* and ^B^). If in Equation (5.45), 
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the terms involving coupling betveen the* general! sed coord in )tco 

2 

are neglected (that is, terms involving S ± B 1f B^ B^ and 
tion (5.13). 


B^), then, the resulting equation has the form of Equa 


If the above analysis stcirting with Equation (5.36) 
is carried out for the case of f.lat plates, the modal equation 
analagous to Equation (5*43 ) which is valid for thin elastic 
shells can be shown to be:- 


{dip B; +(l2)ip + «s)ip Bi + SjOAi, £ Bj +S j E k a 8 X )kl( (B i B J 8 k )]=0 


44) 


where, p as before takes on (one at a time) each one of the 
values that i takes and ,1 £ k ■} i. Equation (5.44) represents 
as many equations as there are p's or i's. Comparing Equa- 
tion (5.44) with Equation (5.43), it is seen that in the case 

of flat plates, the nonlinear te:.*ins of the quadratic type (that 

2 

is, terns proportional to B i and B . ) are absent and only non- 

* 3 2 

linear terms of the cubic type (that is, b£, Bj and B^ B^ B^) 


are present. A similar observation was made earlier with 
reference to Equations (5.13) &n-l (5.20). 


A different approach to obtain the modal equations in 
possible. This involves writing the expressions for the kinetic 
energy T and the strain energy U of the elastic body. The 
expression for U can be obtained from Equations (2.6 )» (2.7) 
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and (2,9) wherein the appropriate relations for the median 
surface strains and CO and for the curvatures 

22 113:0 UFied * ^ iU3 •• u can be obtained in 


terms of the displacement components u, v and w. Assuming 
for u, v and w the product type of functions involving the 
generalised coordinates jj3ee for example, Equation (5.1) ] 
and using the Lagrange’s equation (5.25 ) f the modal equations 
involving the generalised coordinates can be determined. In 
« this approach, the compatibility equation (4.12) and the sim- 
plifications and approximations embodied in it or in the 
equilibrium equation as given by Equation ( 4 . 13 ) are not 
involved, However, in this approach, a set of coupled non- 
linear ordinary differential equations involving the genera- 
lised coordinates will have to bo solved. Th;\s approach of 
deriving the modal equations is similar to what Haters and 
Y/renn [<j] employed in the solution to the problem of non- 
linear f.exural vibration of thin circular cylindrical shells. 
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i) 


( 


( 


jS oct.i on 6. 


In this section the problem of buckling (linear 

> 

problem) and postbuckling of plates and shells arc discussed 
from a general point of view* In Section 4, the governing 
equations for the pootbuckling problem of thin elastic chslls 
have been obtained and they are 1;hc equilibrium equation (4,11) 
and the compatibility equation (4.7) or (4.12). It is helpful 
to reproduce these equations hero. 


D(v 4 W)-(^i + -N ll X ll -N 2J » 21 -2N li « l2 = Z 


and, 


V‘F „ Eh{-tH«) , (^Xl?)} 

In Equation (6.1), — (— + = i) (F) 

\ r i / 


( 6 . 1 ) 


( 6 . 2 ) 


Y or the linear or classical buckling problem, the 
values of the stress resultants ^ and N just, prior 
to buckling do not depend upon tne buckling node V7 so that 
they are determined from the linearised version of the com- 
patibility equation (6.2) wherein the torn b(W ) which depends 
upon the buckling mode XL, is neglected. Thus the stress 
resultants for the lineur budding problem arc determined 
from the compatibility equation « 0. For the pool - 

buckling problem, the stress reoiltanto appearing in Equation 
(6.1) are determined from the compatibility equation (6.T ) 
which is nonlinear. 
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The governing equation* for flat plates can be 
obtained froxa Equations (6.1) erd (6,2) by recalling that 
for flat plates, the principal radii of curvatures ^ ane 
are both ixifinity so that the operator I) vanishes, as 
can be seen from Equation (4,6). 

As in the last section, the function \! is assured 

as a series so that 


W(«,,s 2 ) =£< B 4 \\ s 2 ) (6.3) 

where B^s are coefficients to te determined and W^s arc 

functions each of v/Jiich io assumed to satisfy all the bo\m- 
dary conditions (both geometric and natural) of the problem. 
With the assumed raode shape given by Equation (6.3), the 
compatibility equation (6.^) is satisfied if F^,^) 
is assumed ass 


F(*i.=i) * {t , *F < ,(^,« 1 )+S,.[B 4 F < (S l ,5 1 ) +pf fJ+1.5. B.Bj F*F*} 


(6.4 ) 




where, p* is the applied load (other than any lateral load Z) 

* 

that gives rise to the buckling problem. When p « 0, 
Equation ( 6 . 4 ) reduces to Equation (5.29 )• The set of func~ 
tions {FjJ , {f’-J and arc determined from Equations 

(5.38), (5*39) and (5.40). The function ? 0 appearing in 
Equation ( 6 . 4 ) is determined from the biharaonie equation 




V 4 F„ =0 


( 6 . 5 ) 
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For the lincfir buckling problem, it has been mentioned th;rt 

I 

tho compatibility equation is v F “0 so that, Equation 
(6.4) simplifies to 

F|„.« * f*FoU„s t ) (C.6) 

' v/hcre, p* is the critical value of the load at which buckling 

* 

occurs so that p is the quantity to be determined in the 
linear buckling problem* 


The stress resultants N^ f and If are deiined 

in terras of the stress function F by Equations (4*9) flora 
which we write (in view of ds^ * A^ da^ and dn^ « A 0 drv)s 



r aV i /3 a,\ af i ^ 3 A 2 n af| (6.7) 

^ -35) 3sj Aj \dSi / 3S| A2^3$i / 352- 


The curvature expressions given by Equations (4»4) can be 
written as: 



( 6 . 6 ) 


i — I i I 
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If the same approximation as used in the prcviovc 

sections, namely, neglecting terms involving products of the 

dA« BA . 

gradients of the Lame* paramctci 8 ouch ao (-- - i ) and (-*-'•), 

CCg OD-j 

is made use of, the stress resultants and curvaturec as given 
above can be simplified to: 


and, 


*• *-(£) !.*-■ - : «.-sa 


( 6 . 9 ) 


( 6 . 10 ) 


Using the above relations, the equilibrluo equation 
(6.1) can be v<rltten as: 

, . o .a g • i S ©- S (g) 


Employing, as in the previous sections, the Galerkin *s method, 

we get, 

X D-(W,P)] Sw dA = 0 

A 

where, dA *= ds^ dc 2 is the elemental area on the shell 
median surface. 


With 6w » (6B p ) Wp , . ,« 


' ;he last equation yieldo 


^ t- (w, f)] wp dA = o 


( 6 . 12 ) 
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By meals of Equation (6.3), v/e can v.Tite from Equation (6.11) 


U (W, f ) = [i o + (1,\ 0i + (Vi Q* ♦ ( I,* B? + Z, ( !„>, B, B, + E, % B? 6j 


♦5 *%*} = 0 

(6.13) 

where $ * k 4 i and I Q , (1^)^ etc. are defined as follov/o: 



^ijk” ’ tasf^ F J } 


^ ( r * r») **** 

as? asy* ' dSf 


-2 


dS ( 3Sj 


(F* 



1 

as, as* ]• 


(f».w) 


Define, 


v* 4 i» w i» <ja 

U t \ f = 5 a (I,\ Wp dA 

Wt\ P * i. (I,) < W'pdA. 

p " ( 13)^ Wp d A 

^ ^k\i dA 

/ 

\ 

^A^ijk^> dA 

• 

Note that j f. k f i 





-* 
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Using Equations (0.13), (0.14) and (0.13) in Equation (0.1?), 
we finally get: 




•.E j (y ijf B?6 j+ £ J I k fj c ) jjk /6 l S J B l )] = 0 


(6.16) 


Since p can take on (one at a time) each one of 
the values i can take on, Eqi’ation (6.16) represents t. a 
many equations as there are p's; that is, i'n. This 6tt 
of algebraic equations is nonlinear becaucc each equnticn 

7 p 

has cubic terms (B?, B . B. and B. B 3 terms) and also 

1 1 J i 1 k 

2 

quadratic terms (B^^ and B.^ B^. terms). Furthermore, each 
equation has a right hand side given by J Q ^. J vanishes 
for a flat plate (r^ •* « ard r^ - 00 ) provided the lateral 


loading function Z * 0; in such a case, the set of algebraic 
equations has no constant term. It is to be noted that when 

*s O, Equation (6.16) is similar in form to Equation (5.35) 

• • 

wherein the inertia terra is dropped and B ^ is treated 
as a parameter independant of tine. For the post buckling 
problem of the shell, we have t:> solve simultaneously a set. 
of cubic equations in the coefficients B^s given by Equation 
(6.16). These cubic equations lave a constant term even 
when the lateral load Z = 0, a 3 can be seen from the defini- 
tion of J op given by Equation (6.13). 
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For the linear buckling problem of the shell where 
the buckling loads p* is to be cetermined, Equation (6.16) 
is linearised to give 


ip D i + j o P = 0 (6.17) 

The above represents a set of linear algebraic equations .In 
B^s obtained by \vriting an equat '.on for each value of p. 

For determining the buckling loads p*', the characteristic 
equation is determined from the condition that the coeffi- 
cients of the assumed buckling mode shape become arbi- 
trarily’ large. Thus, the condition becomes: 



( 6 . 10 ) 


which or. expansion gives the following determinantal 
equation: 


<2$ 






0 * ? 


1>11 


(J 1> 31 

( V 41 ... 

(J 1 ^ N 1 


i \z . 

^ j 1 \z 

(J 1 >52 

(J '+2 - 

(J i ha 

= 0 

1 \y 

(J 1 k • 

( v 33 


(J i > N3 

(fa. 19 ) 


• • • • 


It is of interest to sp3cialisc Equation (6.16) 
for the case of flat plates for which r 1 - °® and r - ». 

^ i 


: 01 




C 


Thuc if Z =* 0 so that J 0 p u °» "the governing algebraic 
equation for the post buckling behavior of flat plates is 
obtained from Equation (6.16) a*,: 


b, [a,), p +cy ip b 4 + (j } ) ip af +>- (J 4 ) ijp Bj 




( 6 . 20 ) 


Since in general* B.^ £ 0, the above reduces to a quadratic 
equation, viz., 

{QAp B’ + Zj C\\ JP (B t - B,) ♦r,Z k C^) ykr (B J B k ) 

♦ «4>ip B i + ^ OAjp Bj + CJ,\ p } = 0 (6.2, ) 

Thus, we conclude that for the case of postbucklirg 
behavior of flat plates, the set of algebraic equations in 
♦-he coefficients B^'s to be solved simultaneously consists 
of quadratic equations in B^'s whereas for the case of the 
postbuckling behavior of thin shells, the set of algebraic 
equations to be solved simultaneously consists of cubic 
equations in B^'s and this situation is true when Z = 0. 


It has been mentioned that for the linear buckling 
problem, the stress function determined from the compatib ‘.lity 
equation y^l ? 0 s 0 io given by liquation (6.6). For the post- 
buckling problem, the stress function F is determined from 
the compatibility equation (6.2) which is nonlinear. Bas’d 


on thi? nonlinear compatibility equation and the assumed 
postbuclrling code given by Equation (6.3), the set of 
coupled nonlinear algebraic equations in the coefficient -3 
B^b are given by Equation (6.16), This act consists of 
cubic equations in B^s. An approximation to this set of 
equations can be obtained by working with the linearised 
version of Equation (6.2), namely, 


V 4 F = - Eh D(v/) 


( 6 . 22 ) 


The last equation involves the simplest form of coupling 
between W and F that is possible for the compatibility 
equation. The last equation together with Equation (6.3) 

can be shown to lead to the stress function F given by 

/ 


F(s„s.> ={h*F 0 (s„5 2 ) + £ i B ( F t (s,,s,)} 


( 6 . 23 ) 


It is readily seen that the above is a special case of .? 
given *>y Equation (6,4 )e Using the last equation in the 
equilibrium equation (6.1) and by means of the Galerkin*j 
method, the set of algebraic equations in the coefficients 
B^3 can be shown to be 


J op + ^{< J .><p B i + (J 4p B, ? * Zj (Vijp (°i 8 j>) = 0 


(6.24) 


The above set consists of quadratic algebraic 

equations in the coefficients B^s. In fact, Equation (3.24) 

can be readily obtained from Equation (6.16) if the cubl? 
,32 

terras 'B it 13 ^ and B^B^ terms) arc neglected. 


The 3Ct 
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of equations corresponding to linear buckling pthnt in, 
Equation (6.17)Jcan be obtained from Equation (6.24) if 
the nonlinear terms (B^ and B i terms) are neglected# 
Therefore, Equation (6.24) which is based on the compati- 

4 

bility Equation (6.22) v/hich is linear in \l and F may be 
considered as the "first order postbuckling theory’ of 
shells. Equation (6.24) which is valid for shells can bo 
specialised for the case of flat plates. However, for flat 
plates, it must be remembered that the operator !)(-*) 
vaniehee so that Equation (6.22.1 reduces ; :o y /r F « 0 which 
is the governing compatibility e quation of linear buckling 
theory. Thus, for flat plates, there doc3 not enist any 
"first order postbuckling theory we have only the linear 
buckling theory and the postbuckling theory based on Equa- 
tions (5.1) and (6.2) wherein I»(" " j *= 0. 

If to the governing eqiations for the postbuckling 
problen of thin shells, the effects of damped vibration tire 
included, Equation (6.1 ) ha3 to include the terms pw and Cw 
where p is the mass per unit area of the shell and C is 
the viscous damping coefficient. Thus, Equation (6.1) becomes: 


L(w,F)E [pw +CW +DCV 4 W)~(^- + ) 

L. \ i| •! / 

— N,| 2t|j — N 22 ^22 “ 2 N|2 ^12 “ Z J - 0 


(6.25) 


The compatibility equation (6.2! remains unaltered. In the 
above, Y/ » W(s^,s^,t). Introducing the series for W given 
by Equation (6,3) where B i now is a generalised coordinate 
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r 


and hence a function of timo t, X ( \7 , F ) can be shown to be: 
Lfw,F) = X- (pvv { )E3 { + (ewi) 4- I (W,F) 

(6.26) 

where, I(W,F) is given by Equation (6.15). Thus the modal 
equation obtained by using tho Galer kin's method as in the 
last Ejection becomes: 


ZiKCMW + a\ p 04 + (*>* Bf +y 5 \ p B> 

+f|SMp 04 Bj + <Vo r 3 Bj] (Vijkf 04 Bj B k } ^«P = 0 


( 6 . 27 ) 


In the last equation, J^ |; *» f f M dA and the other 

1 A * 

coefficients such ao (J^)^, iJg)^ 80 cri 8X0 ^^ied 

by Equation (6.15)* If the set of functions { are 03'tho- 

gonal to that f \J \7 dA » 0 for i f p, the coefficients 

A 1 P 


0 for i £ p. It is important to remember that as one 

goes from the large amplitude vibration problem to the post- 
buckling problem, the signs of the stress resultants U , 

N 0 g and JX have to be changed. Therefore, there will be a 

difference in the sign between (F^)^, and (Fj.) where ‘ bhc 

Bubscr j pts pb stand for the postbuckling problem find fv 
stand for the flexural vibratio.i problem discussed in the 




last section. Thus, (7^) ^ Similarly, 


* _* 


<*lV - < F i >tr and (F i typt - - (r r ^>rv* 


# 

If tlie loading functions p and Z are zero, 
and the damping effect is neglected (c =* 0), Equation (£.27) 
reduce:i to the model, equation dsveloped in the last Section, 

namely, Equation (5 • 35 ) • It should be recalled that wh«n 

* 

P * 0 and Z *= 0, the coefficient *= 0. Also, if I ^ is 

•• • 

independant of time t so that *= o and B i » 0, 

Equation (6.27) will reduce to Equation (6.16). Thus, 
Equation (6*27) is a very general equation from which tie 
modal equation for the large amplitude flexural vibraticn 
of shells, plates and beams can be obtained. Also, the 
governing algebraic equations of the postbuckling problem 
of shells and plates as well as of the buckling problem 
(linear problem) of shells and plate? can be obtained fj ora 
Equation (6*27). 


Conclusions 


The general equations governing the large amplitude 
flexural vibration of any thin elastic shell have been derived 
in this report using curvilinear orthogonal coordinates. They 
are two coupled nonlinear partial (ifferential equations and 
they can be readily specialised for shells of revolution or flat 
plates. Although the material is us sidled to be isotropic ioid 
to obey tte Hooke's law for the two dimensional case, tho g«?o- 
mctric norlinearity renders the governing equations nonlinear. 

Based on a series expansion for tho assumed mode 
shape, the modal equations are derived for the case of any thin 
flat plate; and thin elastic shell. It is shown that tho modal 
equation ;for flat plates is of a similar nature as the dynamic 
equation of a mass-spring combination in which the restoring 
force of *;ho nonlinear spring is a cubic and an odd function 
of the displacement. Thus, the nonlinearity associated with 
the large amplitude flexural vibration of thin flat plates of 
any planforra (rectangular, circular, elliptical and so on) is 
of the hardening type and such a conclusion is also true in 
the case of straight beams. On tfco other hand, the modal 
equation obtained for the case of any thin shell indicate p 
that a softening type of nonlincavity is quite possible in the 
case of thin elastic shells. Tho modal equation for the shell 
can be readily specialised to yield the modal equation for a 


shell of revolution or for a flat plate from which the* mod ;1 
equation for straight beams can then be obtained. The contu- 
sions drawn above help to explain the reasons why in the u Girt- 
ing literature, most of the available solutions for the large 
amplitude flexural vibration of thin shells are found to give 
rise to nonlinearity of the softening type whereas in the case 
of flat plates and straight beams, the nonlinearity is always 
found to be of the hardening type. The difference between 
the hardening type of nonlinearity in the case of flat pistes 
and straight beams on the one hand and the softening type of 
nonlinearity in the case of rlng3 and shells on the other hand 
can in general be traced to the curvature (or lacl: of it) of 
the unde.formed median surface of the structure concerned. 

When deriving the governing equations for the lz rgo 
amplitude flexural vibration of flat plates and thin shells, 
the governing equations for the postbuckling behavior of flat 
plates and thin shells can be obtained as a by-product. S udy 
of these equations governing the postbuckling behavior of 
plates and shells leads to the following conclusions. In the 
case of the postbuckling behavior of flat plates, the algebraic 
equations in the coefficients of the assumed series for the 
postbuckling mode are a set of coupled, quadratic equations 
in these coefficients, whereas ii the case of the pootbuc iling 
behavior of thin shells, the corresponding set of algebraic 
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equations are cubic; this situation is true whoa the lateral 
load in r. function Z = 0. If 7> M 0, the set of algebraic 
equations in the case of both plates and shells consists of 
cubic equations in those cocffic Lents. In general, the 
concluo..onc applicable for flat plates also hold true for 
straigh*; beams and those applicable for shells hold good 
for rin^s. 
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